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The spin-gap phase of the one-dimensional t-J-J' model is studied by the level-crossing of the 
singlet and the triplet excitation spectra. The phase boundary obtained between the Tomonaga- 
Luttinger and the spin-gap phases is remarkably consistent with the analytical results at the 
J, J' — > and the low-density limits discussed by Ogata et al. The spin-gap phase has a single 
domain in the phase diagram even if the spin gap opens at half-filling. The phase boundary 
coincides with the K p — 1 line where the Tomonaga-Luttinger liquid behaves as free electrons, 
in the low-density region. The relation between our method and the solution of the two-electron 
problem is also discussed. 
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The discovery of quasi one-dimensional (ID) cuprates 
has stimulated the study of ID strongly correlated elec- 
tron systems. The existence of a spin gap is an important 
factor to be considered in the investigation of supercon- 
ductivity. The effect of frustration was considered to be 
a candidate causing a spin gap in the realistic parameter 
region. However, it was very difficult to determine the 
phase boundary between the Tomonaga-Luttinger (TL) 
and the spin-gap phases, due to the singular behavior of 
the spin gap near the critical point. Recently, the au- 
thor with Nomura and Kitazawa have developed a new 
method to treat the spin gap, that is, the singlet-triplet 
level-crossing method with the twisted boundary con- 
ditions^ In this letter, we will apply this method to 
the ID t-J-J' model and show that the result is consis- 
tent with the exact results in the low-density and the 
J, J' — > limits, and we will clarify the spin-gap phase 
of this model. 

The Hamiltonian of the ID t-J-J' model is written in 
the subspace of no doubly occupied sites as 

S i+ i - riim+i/A), (1) 

1 = 1,2 i 

where = J, j( 2 ' = J'. We also introduce a pa- 
rameter a for the strength of the frustration given by 
a = J' /J. At half-filling (n = 1), this model becomes an 
5 = 1/2 frustrated spin chain. In this case, the ground 
state at a = 1/2 is the two- fold degenerate dimer state 
with a spin-gap, and the ground state energy density is 
— 3/4jJ3'tJL$ Okamoto and Nomura have argued, using 
the level crossing methad, that the fluid-dimer transition 
occurs at a c = 0.2411.0' Upon doping of holes, the sys- 
tem may become metallic, and the spin gap is reduced© 
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but persists for the finite doping. The phase diagram of 
this model for n ^ 1 at a = 1/2, using the exact diagD- 
nalization, was obtained by Ogata, Luchini and RiceLP 
They also made an important achievement regarding the 
critical points for some limits. 

We briefly review their method to obtain the critical 
points in the following two limits: in the limit of J, J' — > 
0, the spin part of this model can be mapped onto the 
case of n = 1, using the factorized wave function,L}EP 

Jeff = J (n.in i+ i) SF + J'(m(l - n i+ i)n i+2 )sF, (2a) 

J' cS = J'(nini +1 ni +2 )sF, (2b) 

where (• • -)sf indicates the expectation value of the non- 
interacting spinless fermion. The effective ratio of the 
frustration a c g is then obtained as 



a e ff(n, a) 



(1 + l/a)n 2 - sf - sl/a 
n 3 - (2sf + sf )n + 2s?s 2 



1 



(3) 



where si = sin(Z7rn)/Z7m. We can obtain the critical den- 
sity n c where the spin gap vanishes, by comparing eq.(|^) 
with tha-jvalue a c — 0.2411 obtained by Okamoto and 
NomuraJf For a = 1/2, n c = 0.7433. On the other 
hand, for the low density limit, the critical value of J/t 
where the singlet pair forms a bound state, can be an- 
alytically obtained by solving the two-electron problem. 
The solution can be derived by adding a Hubbard inter- 
action TLjj = U ^2 i niij7iii to the Hamiltonian, relaxing 
the constraint, and setting U — oo at the end of the 
calculationE' The obtained result is 



■k 
t 



1 + 2a - VI + 4a 2 



(4) 



It is not trivial that this critical value is equivalent to 
the spin-gap phase boundary. We will discuss this point 
later. 

The critical point in the remaining region can be de- 
termined in the following wayfP In this letter we argue 
within the scheme of the bosonization theory. The low 
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energy behavior of lDpfilectron systems is described by 
the sine-Gordon modeMtj 



Ti = Ti. p + TL a 



2fj 



1_L 



dx cos 



a)- (5) 



(27ra) 2 

Here a is a short-distance cutoff, g\± is the backward 
scattering amplitude and for v = p, a 



1 

2^ 



dx 



VyK v 



d9„ 
dx 



dx 



(6) 



where v v and K v are the velocity and the Gaussian cou- 
pling, respectively, for the charge (y — p) and the spin 
[y = a) sectors. In the TL phase (gi± > 0), the pa- 
rameters K a and g\^ are renormalized as K* — 1 and 
gl ± = 0, rcflectinftthe SU(2) symmetry. The phase fields 
are defined adiil-O 



<pv(x),0 v (x) 

in \ - 1 



-ia\p\/2—ipx 



vr{jp) ± vl(p)] Tn Vl m u 



(7) 

\f2lTX 



where v r is the charge (y — p) or the spin (y = a) density 
operator, which satisfy the boundary conditions, 

cj) v {x + L) = <f> v (x) — V^ttUv, (8a) 
6 v (x + L) = e v {x) + V2nm u . (8b) 

These phase fields have the relation [<p v (x), d x 'Q v (x')} — 
iir5(x — x 1 ). The quantum numbers are defined by the 
total number operators (measured with respect to the 
ground state) for right and left going particles (r — 
R,L) of spinaBB 

n* = [{N^ + Nli) ± (N m + N H )}/2, (9a) 
m v = [(N m -N^)±(N Rl -N Ll )]/2. (9b) 

Here the upper and lower sign refer to charge and spin 
degrees of freedoms, respectively. The selection rule 
among them for N = 41 + 2 li: —integer) electrons in 
periodic boundary conditions isHliiP 



m p ±m a 



(_l)n„ 



n„±n a 



(10) 



The finite-size corrections for the excitation energy and 
momentum of the system with length L are described 

by 



E — Eq 
P-Po 



2ttv„ 



2nv n 



L 

2tt 

T 



{s p + s a ) + 2m p k F , 



(11) 
(12) 



where kp = nn/2 is the Fermi wave number, x v = 
{mjK v + n l ?/K u )/2 and s v = m v n v are the scal- 
ing dimension and the conformal spin, respectively, for 
each sector. The corresponding operator is given by 
exp[iv / 2(m.i/0 1 y + n v Qy)\. The boson representation of 
the fermion operator is 



1 



_ e irk F x e i/V2-[r(4> p +crip a )-0 p 



(13) 



where r — i?, L and a =] , [ refer to +, — in that order. 
Next we extract the singlet and the triplet excitation 



spectra [m f 



0,(m CT ,n CT ) = (1,0),(0,±1)]. The 



corresponding operators are \[2 cos \/2<fr a for the singlet 
and \[2 sin ^/2(j) a , exp(+i\/29 a ) for the triplet state.tj' 
Using eqs. (^) and ([l3|), it can be seen that the fermion 
operator should satisfy the anti-periodic boundary con- 
ditions (APBC) for these excitations: ip r ,a(x + L) — 
—ipr !(7 (x). This indicates that APBC plays a pale in the 
elimination of the induced charge excitations.Elf 

The effect of the backward scattering term with the 
coupling gi± is considered using a renormalization group 
analysis under the change of the cutoff a — > e dl a, then 
the scaling dimensions of the singlet (x as ) and the triplet 
(x a t) excitations split logarithmically byJthe marginally 
irrelevant coupling yo(l) = gi±/irv a asEj'EIP 

1 

1 

2 



gi^/irvo- as 
3 yo 



4 y hi L 
1 Vo 



V 



1' 



(14a) 



(14b) 
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where yo is the bare coupling yo = yo{l = 0). In 
the TL region where the backward scattering is repul- 
sive (yo > 0), gi± is renormalized to 0. On the other 
hand, when the backward scattering becomes attractive 
(yo < 0), <7i_i_ is renormalized to — oo and a spin gap ap- 
pears. At the critical point yo = 0, there are no logarith- 
mic corrections. Thus, the critical point is determined 
with high-accuracy by the intersection of these spectra. 
This way of determination of the spin-gap phase bound- 
ary is equivalent to the spin-gap generation in the theory 
of TL liquidsM 

Fig. 1(a) shows the phase diagram of eq.(l) at a = 
1/2. The exponent K p is calculated in the L = 16 system 
in a similar manner as showjaiay-Dgata et al. , which was 
used for the t-J(-J') modelE'LLSEr The phase boundary 
starts from the critical value of the low-density limit, and 
bends at n ~ 2/3. It then flows into the critical point 
of the J, J' — » limit. Thus the spin-gap phase appears 
near half-filling, and in the large J/t region as in the t-J 
model, but it has a single domain in the phase diagram. 
Ogata et al. estimated that the J/t dependence of the 
spin-gap phase in the low-doping region is small, and 
observed the spin gap in the large J./t region at n = 2/3 
using a finite-size scaling methodfl^ These estimations 
are almost consistent with our results. 

Another interesting feature of the phase diagram is 
that the spin-gap phase boundary overlaps the K p = 1 
line where the TL liquid behaves as free electrons in the 
low density region. This point will be clarified later in 
connection with the two-electron problem. 

In the theory of TL liquids, both the singlet and the 
triplet superconducting correlatipns (SS, TS) have the 
same critical exponent 1/K p + l,Ef but TS is dominant 
if the logarithmic corrections are taken into account S3 
However, in the presence of the spin gajx SS is enhanced 
to l/Kp and TS decays exponentiallyLP In the case of 
the t-J model (a = O^pthe spin-gap phase boundary lies 
in the K p > 1 regio: 
hand, at a 

there is no TS region in the latter case 



'■^'(see Fig. 1(c)). On the other 
1/2 it lies in the K p < 1 region, so that 
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Fig. 1. Phase diagram of the ID t-J-J' model at (a) a = 1/2, 
(b) a = 0.2411, (c) a = (TL: TL liquid, SG: spin-gap phase, 
PS: phase-separated state). The spin-gap phase boundary flowp 
into the critical point obtained by the analysis of Ogata et aln 
(plotted by (g in (a)). The contour lines of K p are determined 
by the-d|ata of the L = 16 system, based on the analysis of Ogata 
et aim% 



The phase diagrams for the intermediate values of a 
can be constructed in the same way. As an example, 
we show the phase diagram at a = a c in Fig. 1(b). 
As a is increased from 0, the spin-gap phase boundary 
and the phase-separation boundaryES shift towards the 
small J/t side in the low-density region. On the other 
hand, in the high-density region, the two lines shift to the 
larger J/t side. However, after a exceeds the value a c , 
the spin-gap phase boundary bends and flows into the 
point (J/t, n) = (0, n c ). The phase separation boundary 
also goes to the J/t — ► side following the behavior 
of the spin-gap phase boundary. This behavior indicates 
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Fig. 2. Critical point J c /t versus the strength of the frustration 
a in the L = 16 system at n = 1/2. The critical value is al- 
most constant (J c /t ~ 2.7). This behavior can be explained by 
applying the g-ology. 



that the phase-separated region is always situated on the 
border of the spin-gap phase. 

In spite of the deformation of the phase diagram, the 
critical value J c /t at the quarter-filling (n = 1/2) is al- 
most independent of the strength of the frustration a, 
and is kept at J c /t ~ 2.7 as shown in Fig. 2. Let us 
consider the_reason for this using an argument based on 
the g-ology £3 In order to apply the g-ology, we add the 
on-site Coulomb term Hu to eq.(l) and relax the con- 
straint. The original Hamiltonian is restored when we 
set U = oo. The interaction term in eq.(l) is divided 
into the XY and Ising terms as 

^XY = ^ Y,( S l~ S i+l + S ^ S t + d> ( 15a ) 
i 

W isin g = 2~^2i n it n i+n + nnrii+ii). (15b) 

i 

Since the g-ology is appropriate for the weak coupling 
case, we consider the I = 2 terms of eq. (|l5| ) as corrections 
to the t-J model which belongs to the universality class 
of the TL model. Then their contributions to the g- 
parameters, which are related to the spin-gap generation, 
are identified asEP 

Sgi± = $gi\\ = ~J'(1 + cos4fc F ). (16) 

Note that the ^-parameters are redefined as g^ — g 2 \\ + 
<?2_L — ► 9i\\- For the quarter-filling, eq.(|l6|) vanishes, so 
that the I = 2 terms of eq. (|l5|) do not affect the renor- 
malization flow of the spin part. Thus the frustration 
does not change the critical point at the quarter-filling 
within the scheme of the g-ology. 

Finally, we discuss the equivalence between the level- 
crossing jpethod and the solution of the two-electron 
problem,© and the accordance between the spin-gap 
phase boundary and the K p = 1 line in the low-density 
region. In the low-density limit, the many-body prob- 
lem can be reduced to a two-body problem. The wave 
function of the two-electron system is given by 

* = ^$(*,j)4 T cJj V ac), (17) 

ij 

where $(«,j) = i) for the singlet and — 
— <&(j, i) for the triplet state. It is well-known that for 
a two-body problem the ground state is a singlet state 
as far_as_the bottom of the energy band has no degen- 
eracyO'EjFor periodic boundary conditions, the singlet 
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Fig. 3. Ground state energy (i?o) an d the first singlet and the 
triplet excitation spectra (E B ,Et) of the two-electron system 
with length L = 24 at a = 1/2. i? 3 and Et are calculated us- 
ing anti-periodic boundary conditions. The marked intersections 
will coincide in the L — + oo limit. 



ground state energy is given by Eq = —At + B where B 
is the bound energy, and the critical value where B = 
is given by eq.(|J) without size dependence. However, for 
APBC, the bottom of the energy band is degenerate, so 
that the triplet state can be the ground state. One can 
show that there is no effect of j"> for the triplet state. 
Then the lowest energy of the triplet state for APBC 
is always E t = -At + 8{L) where 6(L) = 8t sm 2 (7r/2L). 
On the other hand, the energy for the singlet state for 
APBC is given by E s = -At + 8{L) + B' where B' is 
the bound energy. Then the level-crossing between the 
two levels takes place at B' = (see Fig. 3). In the 
thermodynamic limit, 8{L — > oo) — and the difference 
between boundary conditions should vanish. This means 
that the critical points for B = and B' = will be the 
same in the thermodynamic limit. In this way, it turns 
out that the singlet-triplet level-cross point is nothing 
but the solution of the two-electron problem (B = 0) in 
the low-density limit. 

In the TL liquid theory electrons behave as free par- 
ticles at K p = 1, so that the spin-gap phase boundary, 
where the bound energy becomes 0, should overlap the 
K p = 1 line in the low-density limit. Moreover, since 
the system in the J < J c region is equivalent to the 
non-interacting spinless fermion system (K p = 0.5), the 
contour lines for 0.5 < K p < 1 will focus on the point 
(J c /i,0). 

In conclusion, we studied the spin-gap phase of the 
ID t-J-J' model using the singlet-triplet level-crossing 
method with the twisted boundary conditions. The ob- 
tained results show remarkable consistency with the ex- 
actly known results. The spin-gap phase appears near 
the half-filling as expected, and also appears as a pre- 
cursor of the phase separation as in the t-J model. The 
spin-gap phase has a single domain in the phase diagram. 
The overlap between the spin-gap phase boundary and 
the K p = 1 line in the low density region is explained in 
the connection with the two-electron problem. 
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